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ABSTRACT 

Through a simple model of particle acceleration and pair creation above the polar caps of rotation- 
powered pulsars, we calculate the height of the pair-formation front (PFF) and the dominant photon 
emission mechanism for the pulsars in the Princeton catalog. We find that for most low- and moderate- 
field pulsars, the height of the pair formation front and the final Lorentz factor of the primary beam is set 
by nonresonant inverse Compton scattering (NRICS), in the Klein-Nishina limit. NRICS is capable of 
creating pairs over a wide range of pulsar parameters without invoking a magnetic field more complicated 
than a centered dipole, although we still require a reduced radius of curvature for most millisecond 
pulsars. For short-period pulsars, the dominant process is curvature radiation, while for extremely 
high-field pulsars, it is resonant inverse Compton scattering (RICS). The dividing point between NRICS 
dominance and curvature dominance is very temperature-dependent; large numbers of pulsars dominated 
by NRICS at a stellar temperature of 10^ K are dominated by curvature at 10^ K. 

Our principle result is a new determination of the theoretical pulsar death line. Proper inclusion of 
ICS allows us to reach the conclusion that all known radio pulsars are consistent with pair creation above 
their polar caps, assuming steady acceleration of a space charge limited particle beam. We identify a 
new region of the P-P diagram where slow pulsars with narrow radiation beams should be found in 
sufficiently large surveys. 

We also show that the injection rate of electrons and positrons into the Crab Nebula inferred from 
the polar cap pair creation model at the present epoch (10^^ electrons and positrons per second) suffices 
to explain the nebular X- and 7-ray emission, the best empirical measure of the instantaneous particle 
loss rate from any pulsar, but that the contemporary injection rate is about a factor of 5 below the rate 
averaged over the Nebula's history required to explain the Nebular radio emission (assuming the Nebular 
radio source is homogeneous). It is not clear whether this discrepancy can be resolved by evolutionary 
effects and by better treatment of nebular inhomogeneity, or is an indication of another particle source 
in the pulsar's magnetosphere. 

Subject headings: Acceleration of particles — pulsars: general 



1. INTRODUCTION 

The region above the polar caps of rotation-powered pul- 
sars has long been the focus of research. From the first, 
e.g. Goldreich & Julian (1969), theoretical papers have 
centered on the idea that a beam of electrons (or ions) is 
accelerated up from the polar cap by the electric poten- 
tial induced by the rotation of the pulsar's intense mag- 
netic field. Shortly afterwards, Sturrock (1971) suggested 
that this particle beam will emit high-energy gamma-rays 
via curvature radiation, which will then interact with the 
background magnetic field and generate copious numbers 
of electron-positron pairs, forming a pair plasma. Rud- 
erman & Sutherland (1975) suggested that this plasma, 
once formed, will act to short out the potential, preventing 
further acceleration. Arons & Scharlemann (1979) then 
showed that this scenario can occur even when the star 
freely emits the particle beam. 
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This basic model has remained qualitatively unchanged 
for the past 25 years. Recently, however, two important 
alterations have been made to the theory. First, the impor- 
tance of general relativistic frame dragging in the creation 
of the accelerating potential was discovered by Muslimov 
& Tsygan (1990, 1992) and later elaborated by Muslimov 
& Harding (1997). Second, several authors (Sturner, et al. 
1995; Luo 1996; Zhang, et al. 1997) have suggested that in- 
verse Compton scattering (ICS) is an important, perhaps 
the dominant, high-energy photon emission mechanism, 
rather than the curvature radiation proposed by the ear- 
lier theories. Zhang & Harding (2000) used a model of a 
Compton-dominatcd pulsar including the GR acceleration 
to obtain predictions for the emitted gamma-ray spectrum 
but did not discuss the pair-production rate itself in any 
detail. 

The dominant emission mechanism and accelerating po- 



Hibschman, 



Astronomy Dept.. 

1 



601 Campbell Hall, U.C., Berkeley 94720-3411; email: 



2 



Hibschman and Arons 



tential combine to set the height of the pair formation front 
(PFF), the point above the polar cap where enough pairs 
have been formed to halt any further acceleration of the 
beam. Given the potential, the PFF height sets the final 
Lorentz factor of the beam, determining the energy budget 
of the polar-cap region. 

This paper calculates, in a direct analj^tic manner, and 
for an arbitrary pulsar, the height of the PFF in the pres- 
ence of the Muslimov and Tsygan accelerating potential 
and including the effects of inverse Compton scattering. 

Qualitatively, two rates combine to position the PFF: 
the creation rate of photons capable of pair-producing 
and the rate at which those photons can pair-produce. 
Section 2 discusses the approximations used for the pair- 
production process, Section 3 discusses the emission rate of 
convertible photons due to various emission mechanisms, 
Section 4 explains the model of the accelerating potential 
used, and Section 5 calculates the resultant PFF height. 

For simplicity, for the rest of this paper we express all 
energies and temperatures in units of the electron energy, 
nieC^ = 511 keV = 5.9 x 10^ K, unless otherwise speci- 
fied. 

2. PAIR PRODUCTION 

The photons capable of pair-producing in the back- 
ground dipole magnetic field are those with energy greater 
than a critical value, 
32 B„ 



- 3 5..1n'A=^^^^^-^'^^^«"^^ W 
where B is the local magnetic field, Eg is the critical mag- 
netic field, Bq = 4.41 x 10^^ G, is the radius of the 
emission point as measured from the center of the mag- 
netic dipole, pe is the radius of curvature of the magnetic 
field at re, and In A is treated as a constant parameter, 
here In A = 20 (Arons & Scharlemann 1979). The nu- 
merical value expresses the magnetic field in units of 10^^ 
Gauss, B = 10^^ Bi2 G, and assumes that that the emis- 
sion point is at = lO^recm and that the radius of 
curvature is a factor fp times the radius of curvature of 
a dipole field on the last closed field line, /?,. = fpPB, 
Pb = 4i?*/36»c = 9.2 X lO'^P^^^cm, where P is the pul- 
sar spin period. The expected possible range of fp is from 
fp = 0.011 P""'^/^, for an extreme multipole field with ra- 
dius of curvature equal to the stellar radius, to fp = 1, for 
a standard dipole field. 

These photons (with e > Ca) must propagate through 
the magnetic field for a distance 

in units of the stellar radius i?* before they pair-produce. 
Since the pair-production opacity is a sharply increas- 
ing function of distance, essentially all such photons pair- 
produce at this distance. This approximation is best for 
low- and mid-range field strengths, where B < Bq/3. At 
higher fields, not only does the opacity increase more grad- 
ually with distance, but other physical processes, such as 
photon splitting (Harding, et al. 1997), may become im- 
portant. Most observed pulsars are in the lower regime of 
magnetic field. 

The created pair will be in a high Landau level and will 
therefore emit synchrotron radiation, bleeding off the com- 



ponent of its momentum perpendicular to the local mag- 
netic field. The synchrotron photons with energies above 
Ca then pair-produce themselves, initiating a cascade. 

The total number of pairs produced by this cascade may 
be found by using a modified version of the generational 
model due to Tademaru (1973). Each generation of pair 
production and subsequent synchrotron emission lowers 
the maximum energy of the photon spectrum by a factor 
of In A (Arons & Scharlemann 1979). Since the expected 
synchrotron number spectrum due to one photon has a 
power-law index of = —3/2, the energy is concentrated 
in the high-frequency end of the spectrum, and we may 
treat all the photons of the i^^ generation as having en- 
ergy ei = In A. If there is no energy loss between suc- 
cessive generations, this produces a final pair multiplicity 
M(eo) = eo/ca from an initial photon of energy cq. 

However, even in the best case, this overestimates the 
pair production, since a significant fraction of the penul- 
timate generation of pairs will have energies less than ta 
and will thereby be lost. Given the power-law index of 
V = —3/2, this reduces the number of photons available to 
the last generation by a factor of Vln A, giving a single- 
photon pair multiplicity of 



M(eo) = 1 + 



1 



(3) 



namely the sum of the single pair generated by the primary 
photon and those generated by the cascade. 

This pair production cascade does not take place instan- 
taneously; it gradually develops over the distance of re/4 
it takes the lowest energy convertible photons to pair pro- 
duce. Most pair production is initiated by photons emit- 
ted close to the surface of the star, so r^/R* ~ 1 and the 
length of the cascade is simply i?*/4. Over this span, the 
pair production rate is linear with distance. Each succes- 
sive generation multiplies the number of generated pairs, 
but it does so over a similarly multiplied propagation dis- 
tance. 

The multiplicity of pairs produced before a height ,s by a 
single photon of energy e emitted at height s,-. then has two 
parts. First, as long as the primary photon has enough dis- 
tance to pair produce, i.e. s — Se > ea/4e, at least one pair 
will be generated. Second, the cascade of secondaries pro- 
duced by that pair will generate e/ (e^Vln A) more pairs, 
reduced by the fraction of the full i?*/4 available to pair- 
produce: 



M(s,; 



-(S - Se) 



(4) 



This assumes that s — Se < 1/4 and ignores the change in 
the magnetic field through the cascade. Since the magnetic 
field decreases with height above the polar cap, this on-the- 
spot approximation slightly over-estimates the amount of 
pair creation, but it does so consistently for all of the emis- 
sion mechanisms considered, and so should not change the 
qualitative results. The quality of these approximations 
will be discussed in a subsequent paper. 

3. EMISSION 

The rate at which these photons are produced depends 
on the emission mechanism. There are three primary emis- 
sion mechanisms which are important over pulsar polar 



Pair Multiplicities and Pulsar Death 



3 



caps: curvature radiation, resonant inverse Compton scat- 
tering (RICS), and nonresonant inverse Compton scat- 
tering (NRICS) in both the Thomson and Klein-Nishina 
regimes. 

Each of these mechanisms produces its own character- 
istic radiation spectrum. However, since all of these spec- 
tra are roughly power-laws harder than N{e) oc un- 
til truncated at some mechanism-dependent maximum en- 
ergy, most of the emitted energy is concentrated in pho- 
tons with energy near this maximum. We then approx- 
imate the spectra by assuming that the entire emitted 
power is converted into photons at the maximum energy, 
so N{e) = Px /ejf , where Px is the power emitted and ex 
the characteristic energy of an emission mechanism X. 

The characteristic maximum energies of these mecha- 
nisms are 

(5) 
(6) 
(7) 



KN 



7 

cr = 2^eB 
3 Ac 3 
2 p 



(8) 



where 7 is the Lorentz factor of the scattering particle, 
cb = B / Bq, p is the local field line radius of curvature, T 
is the neutron-star temperature in units of mc^ , is the 
range of cos^ from which thermal photons arise, and Ac is 
the reduced Compton wavelength, Ac — 3.86 x 10^^^ cm. 
Near the surface of the star, A/i = 1, while at altitudes 
large compared to the size of the emitting region A/x de- 
creases as an inverse-square law. The transition between 
the Thomson and the Klein-Nishina regimes of NRICS is 
at 27TA/U = 1 or 7 = 3.2 X lO^A/x-^Tg-^ Above that 7, 
NRICS is in the KN limit. 

The powers emitted by these mechanisms are 



r>Th 

N 



3.08 X 10^ Tg2A;u(ln(27TA;u) -|- 1) mc^s-\ 



(9) 



27TA^ > 1 (10) 
Pr w 4.92 X 10" 7-iB2^r6ln(7r A/xe^i) mc2s-\ 

JTA^ies' > 1 (11) 

Pc » 5.63 X 10-1^7^8^ rnc^s"^ (12) 
where equations (9) and (10) are from Blumenthal & 
Gould (1970), equation (11) is from Dermer (1990), and 
equation (12) is the standard form for curvature radiation, 
c.f. Jackson (1975). 

Since, for our purposes, only the gross magnitudes of 
these power losses are important, we treat the logarithmic 
terms in equations (11) and (10) as constant for the rest of 
the analysis. Resonant ICS peaks at low energies, so we set 
\n{jTApeg'^) = 1. Non -resonant ICS remains important 
even to high energies, so we set ln(27TA/x) = 4, a value 
representative of the final Lorentz factor of the primary 
beam in pulsars where NRICS sets the PFF. 

For particles moving at moderate Lorentz factors (7 < 

4.6 X IO^Tq^^P^^^), resonant and nonresonant inverse 

Compton scattering are the dominant energy-loss mecha- 
nisms (Sturner 1995), while curvature dominates at higher 
Lorentz factors. For the rest of this section, we adopt a 
model in which the entire surface of the star is hot, giving 
Ap, w 1 for the active low-altitude region. Only NRICS is 



strongly dependent on A/i, and the effects of a hot polar 
cap on that process are discussed in Section 5.1.1. 

At moderate 7, NRICS generally overwhelms RICS. The 
Klein-Nishina NRICS process clearly produces photons 
with the largest possible energy, namely the entire energy 
of the scattering particle. If RICS is to be significant, 
thermal photons must be scattered into resonance with the 
field, so 7T > e_B, therefore we also have ej/^ > cr. In both 
Klein-Nishina and Thomson regimes, then, the NRICS 
photons have a higher energy than the RICS photons, and 
will thus pair produce after propagating a shorter distance. 
If the emission rates of the two processes were equal, the 
higher-energy NRICS photons would set the location of 
the PFF, simply due to their faster conversion into pairs. 

At low 7, the RICS emission rate is greater than NRICS, 
but the RICS rate fades quickly at higher 7. RICS pho- 
tons first pair produce when cr > Ca, or at a minimum 
Lorentz factor of 

if^ «4.8x 104Bf2^pi/2. (13) 

2 €B 



If NRICS produces more photons even at this minimum 
Lorentz factor, it will clearly dominate RICS at all later 

7s. At 7}j™, NRICS operates in the Thomson regime if 

27^'"r < 1, or B12 > 4.02re/^pi/4, and in the Klein- 
Nishina regime otherwise. Most pulsars, then, are in the 
Klein-Nishina regime by the time RICS can generate pairs. 

Dividing the emitting powers by the respective charac- 
teristic energies gives an estimate for the photon produc- 
tion rates for Klein-Nishina NRICS and RICS, 

« 1.54 X 109 7-1t2Ams-i (14) 

iV/{« 1.08 X 10137-2512165-^ (15) 
Given these rates, NRICS in the Klein-Nishina regime 



[tKN 
N 



> 



will produce more convertible photons than RICS {N^ 
Nr) whenever 7 > 7.0 x 10^ B^T^^ A p~\ Clearly, the 
RICS emission is strongest at lower 7. However, even if 
we evaluate the emission rates at 7^*", NRICS produces 
more photons whenever 

Bf^"" <1.9P'/^T^/^Ap^/^ (16) 

For pulsars with magnetic fields less than -Bf^^, RICS 
has no chance of determining the pair formation front; 
at all beam 7's high enough for RICS to generate pairs, 
NRICS has a higher emission rate. 

For fields greater than By/^ , the dominant process is 
still unclear. RICS will only dominate the creation of pair- 
producing photons for a small range of 7, and those pho- 
tons will still take longer to pair-produce than the NRICS 
photons. As the primary beam accelerates through this 
range, RICS will only dominate for a short distance. If 
RICS produces enough pairs in that distance, then it will 
set the PFF, but this is likely true only for high-field pul- 
sars. This expectation is borne out by the more detailed 
calculations in Section 5. 

Curvature emission, on the other hand, produces pho- 
tons with energy large enough to pair produce only when 7 
is on the order of 10^ or higher. At this beam energy, both 
NRICS and RICS emission rates are negligible; the ques- 
tion becomes whether the beam has a chance to accelerate 
to such high Lorentz factors in the presence of radiative 
losses. In Section 5, we find that curvature radiation is 
only important for the short-period pulsars, due to their 
higher voltage. 
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4. ACCELERATION MODEL 

In this paper, we adopt a simple model of the accelerat- 
ing potential, based on the general rclativistic frame drag- 
ging first proposed by Muslimov & Tsygan (1992), and 
similar to low-altitude limits of the potential derived by 
Muslimov & Harding (1997). In this model, the charge dif- 
ference arises because, close to the star, frame dragging re- 
duces the effective rotation frequency fief / to less than the 
observed value. The density of the beam extracted from 
the star equals the Goldreich-Julian charge density at the 
surface, but is insufficient to cancel the charge at higher 
altitudes as frame-dragging weakens and the effective rota- 
tion frequency approaches the asymptotic flat-space value. 
To maintain simplicity, we neglect any additional effects 
of GR, such as those of Gonthier & Harding (1994); these 
effects introduce no more than 10% corrections to our re- 
sults. 

The effective rotation rate is Og// = (1 — KgX~^)fl, 
where x = r / R,, is the dimensionless distance from the cen- 



ter of the star and Kg = 2GI /c^R^ 



0.15/45i?6~' 



with / 



the moment of inertia of the star. The effective Goldreich- 
Julian density is then rjcjix) = TJcji^ ^ KgX^"^), where 
rjQj is the charge density with frame dragging neglected, 
V%j = -n-B/27rc. 

For a more detailed form of the accelerating potential, 
sec Muslimov & Harding (1997). For our purposes, we 
divide the potential into two portions: the low-altitude 
quadratic potential at the surface of the star, representing 
the fringe fields, and a high- altitude linear potential. The 
linear form of the high-altitude potential assumes that the 
distance above the star, s = x — 1, is less than 1; s < 1/2 
suffices. The low-altitude fringe fields occupy s < 6*0 «C 1. 
We assume the cross section of the polar flux tube to be 
round with cross sectional area txO^R^x^ , as is consistent 
with observations of pulsar beaming morphology (Lyne & 
Manchester 1988; Kramer, et al. 1998). 

Due to the conducting boimdary of the surface of the 
polar cone, the high-altitude form of the field is simply 
that of a slightly flared cylindrical charged waveguide op- 
erated below cutoff. Assuming that the beam density 
matches the Goldreich-Julian density at the stellar sur- 
face, the difference in charge density at high altitudes is 
A77(s) = -?,Kgr(i{B{x)/B^)s, whcrc'r]^ = ri%j{R^). Since 
this is above the fringe held zone, the induced electric 
field is almost cylindrically radial, giving a potential of 
$ = nOlRlil - w'^){B^I^'q{s) / B{x)), where 0^ is the colat- 
itude of the last closed field line, 6c = -^/flRt/c, w is the 
fractional position across the polar cap, w = 6^/9c, and 0* 
is the colatitude of the footprint on the stellar surface of 
the current held line. 

Using the potential drop across the polar cap, $cap = 
O^^o, where $0 = is a pole-to-cquator estimate of 

the potential, the potential of a field line at high altitudes 
becomes 

3 

^high{s) = -Kg^capS{l - W^). (17) 

Numerically, $o = 1-26 x 10^'^ Bi2P"^ 1^ and $cap = 
2.63 X Bi2P~^V where wc have neglected the cos a 
inclination term as it is of order 1. This produces an order 
of magnitude stronger acceleration for P ~ 1 sec than the 
previous models of Arons h Scharlemann (1979), which 
depended on the curvature of the field lines to generate 



the starvation electric field. For P < 10 ms, however, 
the frame dragging and field curvature acceleration mech- 
anisms are comparable. 

The low-altitude field can be treated as that of a plane- 
parallel charged slab, so wc expect the potential to increase 
quadratically. The low-level charge density is slightly dif- 
ferent from the Goldreich-Julian density. Following the 
method of Arons & Scharlemann (1979), we choose this 
density difference and the height of the transition between 
the low and high so that the potential and electric field 
are continuous. 

This gives a low-altitude potential of 

^low = ^$oSs'(^^^c(l - w'^f'^ - s) (18) 

where the transition between the low and high forms of 
the potential occurs at s = a/3/2 ^c(l — lu^)-^/^. 

To reduce this to an even more convenient form, we 
drop the cubic term from the low-altitude potential and 

simply extend the quadratic component until it meets the 
high-altitude linear potential. This occurs at 



Sl 



= iy|ee(l - w^f'^ = 8.87 X 10-3p-V2(i _ ^2)1/2. 

(19) 

If we express the particle Lorentz factor in terms of the 
scaled height, t = s/si, we obtain 

Jlowit) = 7it2, ^ ^ I 
lhigh{t) = 71*1 i > 1 

where 



(20) 



3 /3 



2\3/2 



71 = -y -KgtJ-—^ = 5.14 X 104Bi2P-5/2(i _ ^2) 
o V ^ TfhC 

(21) 

Although we have dropped the smooth transition in elec- 
tric field from low to high as an approximation, we are still 
using that underlying model to set the coefficient. The ap- 
proximation is valid as long as the acceleration is not ra- 
diation reaction limited for particle energies « 71 , so that 
7 oc as is true of all of our models. 

This acceleration is strong enough that we can ignore 
the effects of radiative losses on the beam particles due 
to ICS, so this discontinuous electric field causes no dif- 
ficulties. Since the power produced by ICS declines with 
increasing Lorentz factor, once the beam has accelerated 
a short distance above the stellar surface, the ICS losses 
become negligible. Curvature radiation becomes very ef- 
ficient at high Lorentz factor, so if the beam reaches a 
Lorentz factor of order a few times 10^, curvature radia- 
tion will hold it there. We will return to this issue when 
we discuss curvature radiation in more depth. 

For the remainder of the paper, we simply drop the de- 
pendence on w from the potential, in effect examining the 
activity on a typical field line, where neither w nor 1 — w 
is small. 

5. PAIR FORMATION FRONT 

The primary beam from the star will accelerate due to 

the starvation electric field until the beam has produced 
a pair plasma dense enough to short out the field. Sev- 
eral authors have loc;ated this pair formation front (PFF) 
at the point where photons emitted by the primary beam 
first reach an optical depth of one, with regards to the 
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pair-creation opacity. This heuristic was first suggested by 
Ruderman & Sutherland (1975) for curvature-dominated 
cascades, and in that case, it is valid. For curvature radi- 
ation, further acceleration of particles only increases both 
the energy and number of emitted photons, resulting in a 
sharp front of pair production. 

However, if inverse Compton radiation is the primary 
source of convertible photons, this approximation breaks 
down. As the primary beam accelerates, the energy of the 
ICS photons increases, but their number drops. Rather 
than the copious cascades produced by curvature radia- 
tion, where the PFF is a sharp transition, ICS tends to 
produce more gradual cascades where a small number of 
high-energy particles set off a lengthy cascade of pair pro- 
duction. For inverse Compton processes, then, it is not 
guaranteed that the potential will be shorted out as soon 
as the first pair is produced. 

To properly calculate the location of the pair formation 
front in the face of relatively sparse photon production, 
we have to track not only the height at which photons of 
a certain energy are emitted and the distance those pho- 
tons travel before pair producing, but also the number 
of such photons produced and the number of secondary 
pairs produced by the resulting synchrotron cascade. In 
short, we have to integrate the pair production rate until 
we reach some threshold. Since inertial frame dragging 
induces a fractional difference of the beam density from 
the Goldreich- .Julian density of approximately Kg, we set 
the location of the pair formation front by equating the 
number of pairs created per primary particle to Kg. 

5.1. PFF Position 

As discussed in Section 3, for each emission mechanism 
we approximate the emitted number spectrum by assum- 
ing that all of the power goes into creating photons at the 
respective characteristic maximum energy: 



Nx{e) 



ex (7) 



(5(e- ex(7)) 



(22) 



where X E {C, NR. R} , and Px{l) and ex (7) arc again 
the total power and the maximum energy emitted by par- 
ticles with Lorentz factor 7 due to curvature (C), nonres- 
onant ICS (NR), or resonant ICS (R) emission. 

Using the cascade multiplicity (4), the total number of 
pairs produced by one primary particle before a height s 
may be divided into a portion due to the primary photons 
themselves, N±\, and a portion due to the subsequent cas- 
cade of secondaries, N±2, where 



Px{s') 



N±i{s)= r^^" ds'R, . 

cex{s') 

4 f^'"" 
N^,{s) = -= / 
VlnA 



(23) 
(24) 



ds' R,^^^^{s - s') (25) 



where s„ij„ and Smax arc the first and last heights at which 
the primary beam emits photons capable of pair produc- 
ing at or below s. We neglect the variation of with s' , 
making an on-the-spot approximation. 

These heights are the values of Se which satisfy 



with As from equation (2). For fast acceleration, Smin s 
and s — Smax <C s, so these limits are approximately given 
by 

ex(smm) ~ -T- (27) 



~ 12. 

~ 4s 

1 e„ 



s — 



Aex{s) 



(28) 



The lower height Smin is the distance required to accel- 
erate particles to energies capable of pair producing in a 
distance s, while Smax is s reduced by the distance it would 
take a photon emitted at s to pair produce. By the time a 
pulsar beam has reached an altitude where the PFF forms, 
it has typically reached a high energy, so Smax ~ s. The 
form of Smin, however, will depend on the regime of the po- 
tential in which it falls. Since the potential switches from 
its low-altitude form to its high- altitude form at s = s\, 
it is easiest to calculate in terms of the scaled distance, 
t = s/si. 

5.1.1. Non-resonant ICS 

As the beam pulled from the stellar surface accelerates 
to high 7, it will pass from the Thomson regime of NRICS 
into the Klein-Nishina regime. Since this acceleration is 
extremely strong, this transition will occur close to the 
stellar surface, and we can simply use the Klein-Nishina 
form of NRICS for these calculations. For NRICS in the 
KN regime, Pnr{i) is only logarithmically dependent on 
7, so we treat it as effectively constant. Since the power 
is then independent of the Lorentz factor, the NRICS cas- 
cade is insensitive to the details of the accelerating poten- 
tial. 

Furthermore, since the NRICS process generates com- 
paratively few photons, it can only achieve high pair mul- 
tiplicities through extended cascades of secondaries. Using 
the NRICS power, equation (10), in the secondary portion 
of the multiplicity, equation (23), produces a total multi- 
plicity as a function of height of 

Pnr sfRf: 



N±NR{t) = 2- 



CaVlnA 

The critical energy is cnr = 7, 

,/l/2 



-(t-t ■ r 



(29) 



so if we define ,„ = 



tmir, if tm.in. 1 aud tr, 



f' 



ea/4si7it, then tmh 
otherwise. 

Typically, however, the PFF only forms at t » imm, so 
we can neglect t„iin above. The distance to achieve a pair 
multiplicity of n is then 

/ \ 

/ AceaVlnA c \ . , 

The PFF takes place where k = Kg, or 

tPFF,NR = 13.4 B^^'^P^I^T^^fl/^ (31) 

assuming, as before, Kg = 0.15. In terms of the stellar 
radius, this is 



SPFF,NR = Q.ll'd B^^'^P^'^T^^fl 



1/2 



As(ex(se)) 



S- Se 



(26) 



(32) 

For most pulsars, this is well above the regime of the 

fringe fields, t k, \. Likewise, it is still within the linear re- 
gion of the potential (s < 1), provided that B > 2 x lO^'^ 
G. Very young and energetic pulsars may have tpFF < 1) 
as discussed in Section 6. The above formalism remains 
valid regardless. 
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If only the polar cap is hot, the flux of thermal X-rays 
is roughly halved by t = 1 , reducing the number of pairs 
produced. If we model a hot polar cap by applying the 
same technique as above, but attenuating the ICS source 
at t > 1 by l/t^, we find that the pair formation front 
occurs at 



.cap 

^PFF,NR 



or 



41?F,iVfl = 0.399 Sr/PTg-Vp 



(33) 



(34) 



As expected, the PFF occurs at higher altitudes if only 
the polar cap is hot. 

5.1.2. Resonant ICS 

In general, resonant inverse Compton scattering de- 
pends on the large resonant cross section to compete with 

nonrcsonant scattering. The RICS photons arc typically 
lower energy than the NRICS photons and produce less of 
a cascade. As such, the pair production of the primaries 
dominates that of the secondaries, and we use equation 
(24) to compute the multiplicity. 

The steps of computing the PFF height arc the same 
as for NRICS, using the primary-only integral (24) rather 
than equation (25) and t^^„ = ea/8e_BSi7ii. 

Since the emission rate depends on particle 7, the two 
forms of the accelerating potential yield different results. 
There are three major regimes: where the pair formation 
front forms in the low-altitude region, with t,tmin < 1) 
where the PFF forms in the upper region due to photons 
emitted in the lower, with tmin < 1 but t > 1, and where 
the entire process takes place in the high-altitude poten- 
tial, with t,tmin > 1- 

In practice, RICS is only important when it can create 
pairs in the low-altitude potential with tmin < 1. If that is 
so, the number of pairs produced as a function of distance 
in both cases is, to leading order in t„iim 



N±R{t) = 



1 Pk(7i) siR, 



t: 



3 27ies c 
which gives a PFF height of 

tpFF,R = 1350B-//'pV2T6-2/Vp- 



or 



SPFF,R = 12.0 5^2'''' '^6 '^^ \fi 



5.1.3. Curvature Radiation 



(35) 



(36) 
(37) 



Curvature radiation differs from both ICS processes in 
that, as the beam Lorentz factor increases, both the emis- 
sion rate and the average photon energy increase. This 
increase creates a sharp transition in the pair production 
rate once the particle 7 reaches the minimum required to 
pair-produce. Due to this quick generation of a copious 
pair plasma, both the A^±i (24) and N±2 (25) forms of 
the multiplicity lead to nigh-identical results for the PFF 
height. In practice, we find that the approximation de- 
rived from better fits the numerical solution, so we 
will concentrate on that case. 

The curvature power (12) and critical energy (8) give 

for t, tjYiin ^ 1; 



-\jlow 

^^±c — 



2 Pc (71) slR* 



it'°-mtl,„ + 9tZn) (38) 



and for t, t„ 

rhigh 



N. 



in '-^ 

2 Pc (71) slR, 



±c 



(t*5 - 6ttl 



htZ 



(39) 



15e„VhrA 

The mixed case is unimportant, due to the strong t- 
dependence of the multiplicity. 

As before, tmin is the first point at which the beam has 
reached a sufficient Lorentz factor for the curvature pho- 
tons to pair-produce by t, 

Since the power produced by curvature radiation in- 
creases strongly with increasing 7. the pair multiplicity 
is dominated by the emission from the upper end of the 
range, so we can obtain a simple approximation by as- 
suming tmin =0. As most pulsars are in the tppF > 1 
regime, except for the shortest-period millisecond pulsars, 
wc concentrate on high-f, giving 



t: 



high 
PFF.C 



98.2p-2^/^p25/12/l/2 if 



^ PFF.C ^ ^ 



(41) 



When compared to the result of numerically solving for t, 
this is only correct to within a factor of 2-3. This discrep- 
ancy arises because once curvature emission is active, it 
very rapidly shorts out the field, so tppp is close to tmin- 
Neglecting tmin , while giving a straightforward expression 
for the PFF height, introduces an appreciable error. The 
results given in subsequent sections use the full numeri- 
cal solution to the polynomials (38) and (39) for tppp^c', 
equation (41) should be considered a rough estimate. 

In terms of stellar radius, the high-altitude limit be- 
comes 

4i^_^ « 0.678 B-^'^P^^'^^fl'^. (42) 
This depends primarily on the polar cap potential 



as 



SppF CX $cap P 



5/6p-l/12 



Physically, obtaining any pair pro- 



duction from curvature emission requires high potentials 
and beam Lorentz factors. 

These results are inaccurate, however, for the millisec- 
ond pulsars, because radiation reaction prevents the pri- 
mary beam from accelerating to the high Lorentz factors 
needed to generate photons capable of pair-producing in 
the low magnetic fields. The high-altitude accelerating 
electric field is 



^ _ 71 mc^ 



The acceleration due to this field is balanced by the cur- 
vature power loss at a Lorentz factor of 

76a! = 2.26 X 10^ Pjf P- (44) 
Setting the typical curvature photon energy at ^bai to 
^cilbai) = and solving for the magnetic field yields 

P12 = 0.500 P/2/7 (45) 
At fields lower than this, radiation reaction will prevent 
curvature radiation from creating pairs. We cannot yet 
conclude that at higher fields, curvature photons will cre- 
ate pairs, merely that radiation reaction will not prevent 
it. 

If the magnetic field has a radius of curvature of p = P, , 
or fp = 0.011 P~^/2^ then the limiting magnetic field be- 
comes 

P12 =0.137 P*'/'^, (46) 

which would allow more pulsars to pair-produce via cur- 
vature radiation. 
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5.2. Beam Lorentz Factor 
The PFF is located at 

tpFF = T^^n{tpFF,NR, tpFF,R, tpFF,c) 



(47) 



Since the beam stops accelerating at this point, its final 
Lorentz factor is 



Ibeam — Jltppp, iftppF < 1 
"/it PFF, if t PFF > 1 



(48) 
(49) 



Analytically, the estimated final Lorentz factors for each 
emission process, assuming tppF > 1, are 

7jvfl = 6.90 X 106 Bli^p-y^T^'/y^ (50) 

7ii = 6.95 X 10^ B-^'^p-^T-^'^f, (51) 

7c = 5.05 X 106 Bli^p-^'l^-'fll^. (52) 

Since the smallest tppp sets the gap height, the dominant 
mechanism is the one with the lowest predicted beam 7. 

5.3. Polar Cap Heating 

Basic electrodynamics requires that a fraction of the 
positrons created above the PFF be reversed and accel- 
erated back down onto the surface of the star, heating it. 
In the curvature-starvation acceleration model of Arons 
& Scharlemann (1979), this flux is d'^ncjc, while in the 
framc-dragging-starvation model, wc expect a reversed 
flux on the order of KgOcTiGjc, where ncj = rjQj/e is the 
fiducial Goldrcich- Julian number density, see Harding & 
Muslimov (1998); Zhang & Harding (2000); Zhang, et al. 
(2000). 

Since these reversed particles are accelerated through 
the same potential as the primary beam, they strike the 
polar cap with a Lorentz factor ^pff- The downward en- 
ergy flux is then 



' = ^PFprnc frevncjC 



(53) 



where frev is the fraction of the Goldreich- Julian flux that 
is reversed. The energy is then radiated away by a surface 
black body, yielding a polar cap temperature 



^ ^PFFmC^ freyriGjC 



crSB 



1/4 



(54) 



:3.9x 106Bjfp-3/4 (^f^, 



'jPFF'm-c 



1/4 



Since the value of jpff for the ICS processes itself de- 
pends on the stellar temperature, we can solve equation 
(54) for the self-consistent cap temperature and, from 
that, derive the PFF height due to each process. Setting 
frev = Kg^cj this yields cap temperatiircs of 



T6,jvij = 0.730 



rl/6 



n,R = i.38B;,'/''p-'/y'^/'^ 

T6,c = 0.758 <''p-23/48_^i/8 



(56) 
(57) 
(58) 

where we have assumed that tppp > 1, using equation 
(36) for RIGS, the attenuated NRIGS expression, equa- 
tion (33), and equation (49) for the resultant Lorentz fac- 
tor. If tpFF < 1, the appropriate corresponding equations 
would be (36), (31), and (48). For curvature, only the 



high-altitude form is important, equation (41). The low- 
altitude expressions yield 

T^f;^^ = 0.919 <'p-5/i2 / 1/6 (59) 

Tl°^ = 5.13 /3/8. (60) 

The low-altitude forms are most appropriate for the 
youngest, highest potential pulsars. 

These temperatures correspond to PFF heights of 

SPFF,NR,pc = 0.749 f2/3 (g^) 

SPFF,R,pe = 9MB;,''^'p'/^f/' 



(62) 

SpFF,NR,pc = 0.129 P-'^'^'/' (63) 

s'^FF,R,pc = 4.92 B^iP'/' f/' (64) 
Since the curvature PFF height does not depend on the 
temperature, it remains as given in equation (42). 

The regimes in which these expressions apply are some- 
what c:oiiipl(x. The predicted PFF height depends on 
whether the pulsar is in the i < 1 or i > 1 regime. For 
each of RIGS and NRIGS, if t^^^p = s^^^p/si < 1, then 
the low-altitude variant is appropriate, otherwise the high- 
altitude. The shortest PFF height then determines the 
dominant mechanism and the temperature of the polar 
cap. 

However, if the stellar temperature, as derived from a 
cooling model, is higher than this predicted temperature, 
then the polar heating is irrelevant and the normal whole- 
star PFF heights apply instead, equations (32) and (37). 

5.4. Final Pair Multiplicity 

Although the beam acceleration stops at tppp, the pri- 
mary particles will continue to produce pairs even after the 
PFF. After the PFF, cascades fall into two general cate- 
gories: opacity-bounded and energy-bounded. Opacity- 
bounded cascades are primarily limited by the decline of 
the magnetic field with distance from the star; to first 
order, the beam Lorentz factor remains constant. Energy- 
bounded cascades, by comparison, are limited by the 
amount of energy in the beam capable of being converted 
into pairs. Here the beam loses energy faster than the 
magnetic field drops. 

Most pulsars dominated by inverse Gompton processes 
are opacity bounded. By the time the beam has reached 
the PFF, it has typically accelerated to a Lorentz fac- 
tor large enough that further Gompton losses are negli- 
gible. Pulsars dominated by curvature emission, however, 
are energy-bounded. Once the beam has a Lorentz factor 
high enough to make pairs via curvature radiation, radia- 
tive losses quickly slow the beam once the accelerating 
potential has been shorted out. 

For NRIGS from a hot star, the final number of pairs 
generated may then be approximated by counting the 
number of pairs produced by a constant-7 beam, 



X-NR = 



/In A 

Since eais) = elil + s)-y^ 

2 Pnr R* 

knr = 7 - 




R* P< 



NR 



c ea{s)' 
this becomes 



2.1Bi2F-i/2TeV: 



(65) 



(66) 



5eOVlnA 

If only the polar cap is hot, this is reduced by dilution 



to 



cap _ 
^NR 



= 2- 



P, 



NR. 



= 0.094 B12P' 



(67) 



8 



Hibschman and Arons 



using the same approximate attenuation model (cx t ^ for 
t > 1). 

For RIGS, the total number of pairs generated is dom- 
inated by photons which are emitted below the PFF, but 
which arc absorbed above it. The minimum height for pair 
production is t^*" = (€0/27165)", where, as before, a = 1 
for t^'" > 1 and a = 1/2 otherwise. 

The number of pairs produced by the primary photons 
is comparable to that produced by the cascade, so the 
multiplicity is 



If f^'" < 1, the total pair production is dominated by the 




region below t = 1, so t^''^ = 1; otherwise, the pair pro- 
duction is truncated by the decline of the magnetic field 
on the scale of i?* and » s^^. 

This gives two expressions for the final multiplicity, the 
low-altitude estimate, appropriate when B12 ^ P, 

= 7.5 X 10-4 Bli^Tef-^/'^ (69) 
and the high-altitude estimate, appropriate when B12 < 



0.0013 B^aP^/^Tg/, 



1 + 



4.8 + 31nBi2 - 2.5 In P 



VlnA 

(70) 

If curvature sets the PFF, the cascade is energy- 
bounded, and the total number of pairs is roughly the 
overshoot energy of the beam, converted into pairs with 
energy Ca ■ If 72?™ is the minimum Lorentz factor required 
to pair produce, the total number of pairs expected is 



where 



7c 



3 Ac' 



1/3 



PFF 

Ic 



IC 



1.51 X WE 



-1/3 



12 



pl/3j2/3^ 



(71) 



(72) 



This expression for kc depends critically on the slight dif- 
ference between the beam Lorentz factor and 7™*". As 
such, this form cannot be further simplified while main- 
taining accuracy. 

5.5. Death Lines 

Pulsar death occurs when pair-production is too weak 
to generate the multiplicity of Kg required to short out the 
accelerating potential. This sets a minimum potential for 
each process, namely 



^deatfc ^ g Q ^ ;^Ql3 p-5/8^1/2 y 
^nfot = 1-1 X 10"p-"/32/3/8 y 

^R,high 



(73) 
(74) 

2.6X 10"p-*3/18/lV27y (75) 

where wc have used the multiplicity equations (67), (69), 
and (70), along with the corresponding heated cap temper- 
atures, equations (56) and (57). If the star itself is hotter 
than the heated caps, the stellar temperature predicted 
by the cooling model used should be substituted into the 
multiplicity equations, rather than the self-heating model. 

The corresponding death line for curvature radiation is 
best computed by simply finding the minimum $ such 



that the primary beam produces photons capable of pair- 
producing by s = 1, after which point the acceleration 
saturates. Including the decline of the magnetic field and 
the increase in the radius of curvature with increasing dis- 
tance from the star gives a minimum potential of 



= 8.95 X 10i3p-V4_^i/2 y 



(76) 

Radiation reaction, however, also sets a lower limit on the 
required field. Converting the magnetic field limit from 
equation (45) into potential units gives 



death 



c = max(8.95 x lO^^ P" V4/i/2^ 1 32 x lO^^ P" V^A) y 

The resulting death lines are shown in terms of the polar 
cap potential in Figure 1, for both the normal dipole field 
and a p = i?, configuration. As low-to-mid P, NRICS sets 
the death lines, although the dip near P = 0.1 is due to 
the effects of curvature radiation. At high P, RIGS sets 
the death line instead. 

Many of the pulsars near P = 1 are between the two 
death lines. This indicates that either non-dipolar compo- 
nents are important in these pulsars, or that the temper- 
ature of these objects is larger than what the self-heating 
model predicts. We consider non-dipolar effects, namely 
an offset dipole, to be the most likely explanation. 

5.6. Particle Injection into the Crab Nebula 

Study of the emission from the plerions surrounding pul- 
sars can provide empirical constraints on the pair mul- 
tiplicity K. The best studied case is the Grab Nebula, 
energized by its central pulsar, which has a rotation pe- 
riod of 33 msec and an inferred magnetic dipole moment 
H — 3.6 X 10"^° cgs. For this pulsar, our results yield 
K « 5 X 10^. The Goldreich-Julian rate of elementary 
charge loss is Nqj = 2^^V/ec = 1.8 x lO^^ s'^ The 
pulsar's geometry is such that the primary polar cap cur- 
rent extracted by the starvation electric field is made of 
electrons. Then the total rate of injection of electrons 
and positrons into the nebula implied by our model is 



= 2kNgj « 1. 



X 1039s-l. 



The nebular X-ray emission from the equatorial torus 
(Hester, et al. 1995) constrains the electron plus positron 
injection rate from the pulsar into the visible equatorial 
region to be approximately (1 — 2) x 10^^ particles s~^ 
(Kennel & Goroniti 1984; Gallant & Arons 1994, J. Arons 
& A. Spitkovsky, in preparation). The flow feeding the 
equatorial region fills a latitudinal sector with full angular 
width of approximately 26 ~ 20°, assuming the fiow from 
the wind termination shock to the visible torus is purely 
radial. If the pair outflow from the pulsar is spherically 
symmetric, our model of the multiplicity implies ~ 3 x 10^^ 
electrons and positrons per second are being fed into the 
equatorial X-ray emitting region, well more than enough 
to explain the X-ray observations. 

In reality, both non radial divergence of the post shock 
fiow in the nebula and the likely decline of the pair fiux 
toward the boundaries of the polar flow (the region of the 
polar cap outflow which feeds the equatorial) both suggest 
that our estimate of the particle flux into the X-ray torus is 
an upper limit. Therfore we conclude that polar cap pair 
creation in the Grab pulsar is adequate to fully explain 
the present day particle injection rate measure by X-ray 
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Fig. 1. — Pulsar death lines for a self-heated polar eap. The solid line is the death line for a dipole magnetic field, while the dotted line 
assumes a p = i?* multipole field. Pulsars, in the dipole case, dominated by NRICS are marked with circles, those dominated by RICS are 
marked with stars, and those dominated by curvature are marked with crosses. 



observations of the Crab Nebula. Other plerions will be 
considered elsewhere. 

By contrast, it is not so clear that polar cap pair cre- 
ation can explain the total particle injection rate into the 
Crab. It has long been known that the radio emission from 
the Nebula requires a particle injection rate, averaged over 
the Nebula's history, of approximately lO^'^s"^ (e-ff- Rees 
& Gunn (1974)), an estimate probably uncertain to a fac- 
tor of 3. Our calculations yield a total injection rate today 
about a factor of 5 below the average rate required to ex- 
plain the radio emission, assuming no Nebular evolution. 

Whether this discrepancy is a serious problem is un- 
clear. The estimate of the average rate of pair injection as- 
sumes the Nebula to be a homogeneous synchrotron source 
- clearly, it is not. The pair creation rate probably was 
higher in the past, when the pulsar was rotating some- 
what more rapidly, so long as the star was cool enough 
for pair creation to operate on the curvature dominated 
branch shown in Figure 5, as is the case today. On the 
other hand, adiabatic losses as the nebula expands might 
increase the required average rate. An additional magne- 
tospheric source of particles, such as an outer gap, may ex- 
ist. These more sophisticated nebular modeling and global 
magnetospheric questions are outside the scope of our in- 
vestigation, and we leave this subject with the remark that 
polar cap pair creation may be adequate to explain the to- 
tal particle injection rate into this well known synchrotron 
source, but the case is not as clear as it is for understanding 
the injection of the particles that create the high energy 
synchrotron emission. The application of our models to 
other pulsar-plerion systems is also a topic for a separate 
investigation. 

6. DISCUSSION 



A broad survey, such as this, best answers broad ques- 
tions about the population. With these results, we 
seek to answer two questions, namely "Which gamma-ray 
emission mechanism dominates in which region of pulsar 
space?" and "How well does this model explain the ob- 
served distribution of pulsars in P-P space?" . 

For any pulsar, the smallest of tpFF,R, tpFF,NR, and 
tpFF,c sets the height of the pair formation front. For a 
whole-star temperature of 10^ K, assuming no special field 
configuration, the pulsars of the Princeton pulsar catalog 
(Taylor, et al. 1993) are categorized as shown in Figure 
2. Figure 3 shows this same dataset in terms of the more 
physical variables of surface magnetic field and polar cap 
voltage. This categorization depends on the assumed tem- 
perature and size of the thermally emitting region, but sev- 
eral basic conclusions can still be drawn. The largest-field 
pulsars tend to be dominated by RICS, the largest- voltage 
by curvature radiation, and the rest by NRICS. 

As can be seen in Figure 3, the important boundaries 
are between curvature and NRICS at high voltage and 
between NRICS and RICS at high magnetic field. Com- 
paring the predicted PFF heights due to RICS and NRICS 
yields 



tpFF,NR 
tpFF,R 



0.0099 5^2'/ V/^Te^'/Vp"'/'- 



(78) 



The botmdary between NRICS and RICS is almost com- 
pletely set by the magnetic field. RICS will dominate only 
for pulsars with very high fields, B12 > 12AP-^/'^'^Tq^^\ 
Similarly, comparing the PFF heights due to curvature 
and NRICS yields 



.high 

= l.^B-^^^P'l^n, 

tpFF,NR 



(79) 
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Fig. 2. — P — P diagram showing which pulsars have the PFF height set by curvature (crosses), NRICS (circles), and RICS (stars) 
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Whenever this ratio is less than one, the PFF is set by cur- 
vature emission; in this regime, RICS is never important. 
Thus, curvature dominates whenever P < 0.22 bI^^Tq 
This boundary is less sensitive to the magnetic field 
than that between RICS and NRICS, but depends more 
strongly on the temperature. 

Perhaps the more interesting question is where this the- 
ory can comfortably generate enough pairs to short out 
the accelerating potential, assmning that pair creation is 
required for pulsar emission. Since many pulsars are domi- 
nated by NRICS, and NRICS depends strongly on the size 
and temperature of the thermally emitting region, both 
the cooling model chosen for the pulsar and the degree of 
polar-cap heating are extremely important. 

The case of a cool (10^ K) star with the self-consistent 
polar cap temperatures derived in Section 5.3 is shown in 
Figure 4. Using the polar cap NRICS multiplicity, equa- 
tion (67), we find that the cascade cannot produce suffi- 
cient pairs to short out the electric field if the predicted 
PFF is above spFF = 0.25. Given this criterion, this model 
succeeds for pulsars with timing ages less than about lO'' 
years, but fails for older pulsars. 

Figure 5 shows the results of adopting a cooling model 
derived from Tsuruta (1998), which includes the effects of 
the magnetic field on thermal conductivity, but no polar 
cap heating. Since heat conduction is easier along mag- 
netic field lines than perpendicular to them, the polar re- 
gions remain hotter than the standard cooling of an un- 
magnetized pulsar, out to ages greater than 10^ years. Full 
2D calculations have only been extended to ages of several 
million years, but the polar regions show no sign of rapid 
cooling up to that point. For our model, we assume that 
more rapid cooling, with a power-law slope of —0.375 be- 
gins at 10^ years. This power law represents thc^ effects 
of internal heating: with no additional heating, the field 
would decay exponentially beyond this point. 

The precise temperature model used is 

logT = 6.65-0. Hog T, T<Wyr , . 

logr = 8.575 -0.375 log T, r > 10^ j/r ^ ^ 

where r is the spindown age of the star. This is a generous 
model, in that it predicts a hot polar cap to late ages, but 
it demonstrates the strong effect of the cooling model on 
the pair-creation process. 

Using the NRICS multiplicity, equation (66), we find 
that the cascade cannot produce sufficient pairs to short 
out the electric field if the predicted PFF is above s pff = 
0.45. The combination of magnetic cooling and internal 
heating allows pulsars to remain active out to timing ages 
of 10^ years, without further modifications. 

The notable exceptions in both of these cases are the 
oldest pulsars and the millisecond pulsars. The millisec- 
ond pulsars, for the most part, require particles to reach 
a Lorentz factor so large that their acceleration is trun- 
cated by radiation-reaction before reaching the threshold 
for pair-production, while the old pulsars simply do not 
provide enough voltage and field strength to pair-produce. 
These pulsars could all be explained by invoking a mag- 
netic field more complicated than a simple star-centered 
dipole, however. Due to the success of the rotating vec- 
tor model for pulsar polarization (Radhakrishnan & Cooke 
1969), we believe an offset dipole to be more plausible than 
simply adding an arbitrary mixture of higher-order multi- 



poles. 

An offset dipole could both increase the surface field 
strength by moving the center of the dipole closer to 
the surface and reduce the effective radius of curvature 
through gravitational bending of photon orbits. Decreas- 
ing the radius of curvature to i?* allows almost all pulsars 
to generate sufficient pairs, as shown in Figure 6. Any in- 
crease in the effective surface field would produces similar 
results, when combined with the effects of gravitational 
bending of the photon orbits (J. Arons, in preparation). 

For a conservative model of a self-heated polar cap, the 
predicted final pair multiplicity is shown in Figure 7. If 
ICS photons create the PFF, the typical multiplicity is 
low, in the range of 1-100 for the bulk of pulsars, much 
lower than the 10'^ commonly assumed. If curvature pho- 
tons stop the acceleration, however, such high multiplici- 
ties can be reached. 

Since the calculation of the expected multiplicity in- 
cludes the decay of the magnetic field with distance, the 
multiplicity gives a good estimate of which pulsars can 
form the PFF, as discussed in Section 5.5. The total mul- 
tiplicity shows a steady decline with timing age, dropping 
to 0.1 near 10^ years, below which the pair creation is too 
sparse to form the PFF. 

In Figure 8, we show the pair multiplicity using the 
Tsuruta temperature model and no cap heating. The in- 
creased stellar temperature keeps the pair multiplicity ex- 
pcc'tcid due; to NRICS between 10 100 out to a few tc;ns of 
million years, while the lack of polar cap heating results in 
curvature dominance and correspondingly high multiplici- 
ties for some of the youngest pulsars. If polar cap heating 
were included, NRICS would truncate beam acceleration 
for most of those pulsars, dropping the multiplicity back 
to the 100 - 1000 range. 

7. CONCLUSIONS 

Even without considering the competition between the 
various emission processes, the most robust result of this 
work is that most pulsars can generate enough pairs to 
short out the starvation electric field without invoking a 
complex dipole structure, if the polar cap temperature 
is over 10^ K. However, without additional heating, the 
surface temperature drops below this value after c. 10^ 
years, making older pulsars problematic. For the hot- 
ter stars, nonresonant inverse Compton scattering easily 
creates pairs in the simple dipole case, unlike the earlier 
models of pure-curvature cascades (Arons & Scharlemann 
1979). Like any cooling neutron star, pulsars should pro- 
duce observable thermal X-rays, with some enhancement 
due to the small, but hot, polar cap. 

Due to the effects of ICS, the expected pair multiplicities 
are lower than those expected from curvature cascades, in 
the range of 1-100, in contrast to the multiplicities of 10^ 
or higher produced by curvature cascades. Despite the in- 
crease of the ICS cross-sections at low energy, the primary 
beam of particles extracted from the pulsar still acceler- 
ates to Lorentz factors of 10^ or higher, due to the lag 
between the emission of pair-producing photons and their 
actual conversion into pairs. 

The improved death line calculations reveal a regime of 
pair-production available at high magnetic fields and large 
periods, first predicted by Arons (1998). At large periods, 
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Fig. 5. — PFF height vs. timing age, using a temperature model derived from the magnetic cooling results of Tsuruta (1998). Pulsars 
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with crosses. Above the dotted line at sppp = 0.45 i?*, NRICS fails to create sufficient pairs to halt acceleration. 
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the polar cap is small and the accelerating electric field 
weak. The small cap prevents NRICS from producing 

enough pairs, if only the cap is hot, but the slow accel- 
eration allows RICS to do so easily. The small polar cap 
also produces a narrow beam, making these objects diffi- 
cult to detect. The recent discovery of PSR J2144-3933 
by Young, et al. (1999) appears to confirm these results, 
as does the more recent theoretical treatment by Zhang, 
et al. (2000), which should encourage further searches in 
this region of parameter space. 

Since each of the different emission processes dominates 
in different regimes of pulsar space, the radio emission 
process must be insensitive to the precise way in which 



the pair plasma was formed. The comparatively low pair- 
production multiplicities found constrain the radio emis- 
sion process by limiting the available particle densities. 

More detailed studies are, of course, required to under- 
stand the expected pair and gamma-ray spectra from these 
stars. For any given pulsar, more precise calculations are 
required to compute the beam 7 and the PFF height as 
functions of distance from the magnetic axis, but the gen- 
eral trends illuminated here should hold. 

Our research on pulsars was supported in part by NSF 
grant AST 9528271 and NASA grant NAG 5-3073, and in 
part by the generosity of California's taxpayers. 
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